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^\ Abstract 
O 

. We present a classical model for bulk-ensemble NMR quantum computation: 

the quantum state of the NMR sample is described by a probability distribu- 
\ tion over the orientations of classical tops, and quantum gates are described 

by classical transition probabilities. All NMR quantum computing experi- 
ments performed so far with three quantum bits can be accounted for in this 
■ classical model. After a few entangling gates, the classical model suffers an 

Q"^ , exponential decrease of the measured signal, whereas there is no correspond- 

\ ing decrease in the quantum description. We suggest that for small numbers 

of quantum bits, the quantum nature of NMR quantum computation lies in 
the ability to avoid an exponential signal decrease. 



(N 



^3 : I. INTRODUCTION 

The original proposals [0,0 for quantum computing using high-temperature, liquid-state 
^ ' nuclear magnetic resonance (NMR) sparked an explosion of interest. The resulting work on 
■ NMR quantum-information processing has attracted a great deal of attention. In addition 
to a long list of publications [lll-pO|], there have been numerous news reports, especially on 
the recent NMR experiments on quantum error correction ^ and quantum teleportation 

i- 

The fundamental information-processing elements in NMR are two-level nuclear spins, 
called quantum bits, or qubits for short, which are bound together in a single molecule. A 
liquid NMR sample contains a macroscopic number of molecules, each of which functions as 
an independent information-processing system. The molecules are initially in thermal equi- 
librium at high temperature; the nuclear spins are only weakly polarized along the direction 
of a strong applied magnetic field. NMR techniques cannot control the quantum states of 
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individual molecules; instead, all the molecules in the sample are manipulated in parallel. 
Moreover, NMR readout techniques are sensitive to the average magnetization of the en- 
tire sample. For these reasons the use of high-temperature, liquid-state NMR techniques 
for quantum- information processing is called bulk- ensemble quantum computation. Clever 
techniques have been devised to map the coherent quantum-mechanical manipulations re- 
quired for quantum computation to this situation in which one neither controls nor makes 
measurements on individual molecules P|,|^,|5| JTT||r3|] . 

This paper deals exclusively with bulk-ensemble quantum- information processing, not 
with proposals (see, for example, |^|) for using nuclear spins in situations, as in solid-state 



systems, where they can be highly polarized. Although the paper is phrased in terms of 
and aimed at NMR quantum computation, the analysis is not specific to NMR and could 
be applied to any realization of bulk-ensemble quantum computation. 

Bulk-ensemble NMR quantum computation is in principle scalable to many qubits 



but it is unclear whether it will ever lead to useful quantum computations involving more 
than a few dozen qubits [|l],^ . Nevertheless, the experiments performed so far, which use only 
a few qubits, are cited as important demonstrations of the principles of quantum-information 
processing and quantum computation. There has, however, been persistent skepticism about 
whether these experiments demonstrate genuine quantum- information processing. Recently 
these doubts have been brought into sharper focus with the realization that all quantum 



states used in present NMR experiments are separable p3[, i.e., unentangled. Entanglement 
is generally thought to be an essential feature of quantum computation p^-p6 . 

In this paper we explore the "quantumness" of bulk-ensemble quantum computation. 
There being no general definition of what it means to be doing genuine quantum-information 
processing, we begin by proposing three different criteria, which are investigated in the 
remainder of the paper. 

• According to the first criterion, an A^-qubit NMR experiment does not demonstrate 
genuine quantum-information processing if all quantum states that occur during the ex- 
periment are separable (unentangled), i.e., are equivalent to classical ensembles of spin- 
ning tops. The motivation for this criterion is that for separable states, the statistics 
of any measurement performed on the NMR sample can be described in the language 
of classical probabilities for spin orientations. With respect to this criterion, NMR ex- 
periments for up to about 15 qubits do not demonstrate genuine quantum-information 



processing ||2^. There is, however, an immediate objection to this criterion: even 
if all states involved in the experiment are classical, the transformations between the 
states might be essentially quantum mechanical, not having any description in classical 
language. 

Our second criterion addresses this objection by stating that an A^-qubit NMR ex- 
periment does not demonstrate genuine quantum-information processing if we can 
construct an overall classical model for the experiment, by which we mean that in 
addition to the experiment's satisfying the first criterion, the transformation T that 
maps the initial quantum state to the final quantum state can be described by classical 
transition probabilities. The transition probabilities are required to be fully specified 
by T; in particular, they must be independent of the initial and final states. We 
formulate an overall classical model below, which shows that, with respect to this 
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second criterion, NMR experiments for up to about 6 qubits do not demonstrate gen- 
uine quantum-information processing. An objection to this second criterion is that an 
overall classical model does not constitute a valid model for a quantum computation, 
because the complexity of the overall transition probabilities might increase exponen- 
tially with problem size; according to this objection, computational complexity is at 
the heart of the question of whether a given computation is quantum mechanical. 

Our third criterion incorporates this objection as follows. It states that an A^-qubit 
NMR experiment does not demonstrate genuine quantum information processing if we 
can construct a gate-by-gate classical model, by which we mean that in addition to 
the first criterion being satisfied, any gate U implemented in the experiment can be 
described by classical transition probabilities, which, as before, are fully specified by 
U and which, apart from trivial contributions for the qubits not involved in the gate, 



are independent of the total number of qubits We devise a gate- by-gate classical 
model below. It shows that all 3-qubit quantum computing experiments performed 
so far do not demonstrate genuine quantum-information processing according to 
this third criterion. 

In the NMR literature, transformations are often represented by diagrams that depict 
spins as arrows evolving on spheres; see, e.g., ||28| or the "effective fields picture" in ||14|| , 



where a diagram representing a controlled-NOT gate is given. This diagram illustrates what 
is new in our approach. The diagram assumes that the sample is split into two subpopu- 
lations, according to whether the control qubit is up or down. Thus it provides a classical 
description of the gate for initial conditions corresponding to computational basis states of 
the control qubit, but not for arbitrary, possibly entangled initial states. Consequently, it 
does not give a classical model of the controlled-NOT gate in the sense of our third criterion. 

We emphasize that our gate-by-gate classical model is not the same as a simulation of 
a quantum algorithm on a classical computer. First, in our classical model, measurement 
statistics are determined by probabilities for classical spin orientations; in a classical sim- 
ulation, where one calculates the quantum amplitudes in, say, the spin-up-spin-down, or 
computational basis, one cannot go from the amplitudes to measurement statistics using 
classical probability logic. This means that our classical model provides a hidden-variable 



description of the experiment |2^, but a classical simulation does not. Second, the computa- 
tional complexity of a classical simulation increases exponentially with the number of gates 
||24]| , whereas the complexity of our transition probabilities for an entangling gate increases 
only linearly with the number of qubits. 

Equally important is a difference between our gate-by-gate classical model and the quan- 
tum description of an NMR experiment. In our classical model, the strength of the mag- 
netization signal decreases exponentially with the number of entangling gates. Though this 
decrease can be staved off for a few entangling steps, it cannot be put off indefinitely. Such 
a decrease is absent from the quantum description. Thus it is straightforward to design 
experiments that cannot be described by our classical model. Two-qubit experiments that 
fall into this class are and probably also [|10|, although in the latter case this is not clear 
from the published data. In order to prove conclusively, however, that a bulk-ensemble 
NMR experiment demonstrates genuine quantum-information processing in the sense of the 
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third criterion given above, one must rule out all gate-by-gate classical models, not just the 
particular model given in this paper. 

The paper is organized as follows. In Sec. |T|, we give a brief description of bulk-ensemble 
NMR quantum computation, and in Sec. |T|, we summarize the proof given in p3| that the 
states used in high-temperature NMR experiments are classical up to about 15 qubits in the 
sense of the first criterion above. In Sec. we define the transition probabilities that are 
used in our overall and gate-by-gate classical models. Used directly in a gate-by-gate classical 
model, these transition probabilities give rise to a signal decrease by a constant factor at 
each entangling multi-qubit operation. Section [V| shows how to modify the gate-by-gate 
model so as to avoid this signal decrease for a small number of entangling steps. In Sec. |V| 
we consider the implications of our results for bulk-ensemble NMR quantum-information 
processing. 



II. NMR QUANTUM-INFORMATION PROCESSING 

All NMR quantum computing experiments performed so far P|-p!^ work according to 
the following principles. The state of the sample is described by a density operator p for the 
spins (qubits) in each molecule. The molecules are prepared in an initial state 

p= (l-e)M + epi , (2.1) 

where M = //2^ is the maximally mixed density operator for N qubits (/ is the unit opera- 
tor) and pi is a density operator, usually chosen to be the projector onto the computational 
basis state |0 . . . 0). The parameter e scales like 

aN , , 

^=^- (2-2) 

Here a = hu/2kT, where u is the average resonant frequency of the active spins in the strong 
magnetic field, determines the polarization of the sample. If pi is a pure state, p is called a 
pseudopure state [0]. Procedures for synthesizing a pseudopure state from an initial thermal 
density operator are described in [^H,|TT],|13| . 

For typical molecules in present magnetic-field strengths, the average resonant frequency 
is roughly u ~ 300 MHz, which at room temperature, gives a ~ 2 x 10~^. It is difficult 
to know the actual value of a and e in the experiments, because the polarization is not 
absolutely calibrated. In evaluating the experiments, we adopt the scaling of Eq. ( |2.2| ), with 
a conservative value of a = 2 x 10~^, which is meant to take into account inefficiencies in the 
experiments, especially loss of polarization involved in the synthesis of a pseudopure state. 

One interpretation of e is that it specifies the fraction of molecules in the sample that 
occupy the desired initial state pi. This interpretation is not unique, and it is not mandated 
or even preferred by quantum mechanics. It is typically promoted to the level of a physical 
fact by advocates of NMR quantum computation, yet it becomes a physical fact only if one 
actually prepares a fraction of the molecules in a particular state ||2^, a situation that does 
not apply to a high-temperature NMR experiment. The freedom to interpret p in terms of 
other ensembles underlies the conclusions of this paper. 

After synthesis of the desired initial state, a sequence of radio-frequency pulses, alter- 
nating with continuous evolution, is applied to the sample. We first consider the case where 
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the evolution is described by a unitary transformation U. Applying the unitary operator U 
to p results in the state 

p°"* = UpU^ = {1- e)M + eUpiU^ = (1 - e)M + ep°"* . (2.3) 

The totally mixed part of the state is unaffected by the unitary transformation. The out- 
put state retains the form with the same value of e, but — and this is the essence of 
the bulk-ensemble paradigm for quantum computation — pi undergoes the desired unitary 
transformation. 

The same conclusions hold for transformations implemented using gradient pulses to- 
gether with diffusion in the sample ; these are equivalent to mixtures of unitary transfor- 
mations, as in 

pOut ^ Y^PiUipUl = (1 - e)M + eY^p,UipiUl = (1 - e)M + ep^ , (2.4) 

where the pi > are probabilities. We assume that decoherence also preserves the form 
of the density operator, leaving e unchanged. This assumption is certainly justified when 
decoherence is simulated by unitary transformations, as in the gradient pulses used to sim- 
ulate decoherence in 0. It is also true for naturally occurring decoherence for times short 
compared to the characteristic time for relaxation to thermal equilibrium. These are the 
times of interest for quantum- information processing. 

The last step in any NMR experiment is the readout. By applying radio-frequency 
pulses and then measuring the transverse magnetization of the sample, one can determine 
all expectation values of the form 

tr(p°^^V^, ^...(^ap^)=e tr (p°"V^, ® ■ ■ ■ ® (x^^) . (2.5) 

The tensor product in this expression includes one operator for each spin; a^^, denotes the 
unit operator / for the kth spin if jSk = 0, and it denotes a Pauli matrix if = 1, 2, or 3. 
In writing Eq. ( p.5|) and all such expectation values in this paper, it is assumed that there 
is at least one Pauli matrix in the tensor product (not all the /5's are zero). The maximally 
mixed part of the density operator does not contribute to the measured expectation values, 
which are determined by the state p°"* that undergoes the desired evolution. The parameter 
e appears naturally as a measure of the strength of the magnetization signal (or of the 
signal-to-noise ratio). 



III. SEPARABILITY OF STATES USED IN NMR 

Our conclusions rest on the freedom to write states of the form ( |2.1| ) in terms of proba- 
bility distributions over spin orientations for N classical tops. We review one such represen- 
tation [^], the foundation for our work, which provides a classical description for all states 
of the form ( p.l| ) provided that 

1 , , 

' - = 1 + 2^^-i • ^^-^^ 

If we assume that a = 2 x 10^^ in Eq. (|2.2|) , this inequality holds for < 16 qubits. 
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The argument is straightforward [|23|. The most general pure product state of qubits 
has the form 



Pi(n) = Pi(ni, ...,nN) = + ni ■ a) ® ■ ■ ■ ® (/ + n^v ■ <?) , (3.2) 

This state can be interpreted as classical tops pointing in the directions given by the unit 
vectors ni, . . . ,nAr, denoted collectively by n. Since the operators Pi{n) form an overcom- 
plete basis in the space of linear operators acting on qubits, any A^-qubit density operator 
p can be expanded as 

p = j dVLnWp{n)Pi{n) , (3.3) 

where dVLn = dVl^^ ■ ■ ■ dQfi^ . The overcompleteness means that the expansion coefficients 
are not unique. One choice is pc 



w'{h) = ti[pQ,{n)) , (3.4) 
where the operators Qi{h) are defined by 



Qi(n) = j^il + 3^1 ■ a) ® ■ ■ • ® (/ + 3niv • ^) ■ (3.5) 



The expansion coefficients wf{n) can be positive or negative, but they obey the bound 

/ smallest \ 



eigenvalue I = ---^ . (3.6) 
V ofgi(n) J w 



When w''{n) is everywhere nonnegative, it can be interpreted as a classical probability 
distribution for the A^ spins to point in the directions ni,...,njv. Since the maximally 
mixed density operator M has probability density w^\n) = 1/{Att)^, it follows that for a 
density operator p of the form ( |2.1| ) with e < rj, 

^'{n) = ^ + {n) > > . (3.7) 

A density operator for a joint system is separable if it can be written as a nonnegative 
linear combination of product density operators. When w''{n) is everywhere nonnegative, the 
expansion ( p.3| ) provides a separable representation for p. Related work on the separability 
of states near the maximally mixed state can be found in p9| and 0]. Separable states 
have no quantum entanglement The importance of separability in this paper is that 

a separable state of A^ qubits can be interpreted in terms of an ensemble of classical tops, 
because the expectation values ( p^.5|) have the standard form for an ensemble with probability 
distribution w^{n): 

ti(^paf3^ (g) ■ ■ • (g) 0-^^) = / dnnWp{n) tr(^Pi{h)af3^ ® ■ ■ • ® cr^^^ . (3.8) 

. ' 

= (ni)/3i ■ ■ • (n7v)/3^ 
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In this expression, {nj)^. = 1 if j3j = 0, and {nj)f3. is a Cartesian component of the vector 
fij if Pj = 1, 2, or 3. 

If one adopts our first criterion, none of the NMR experiments performed to date has 
done any genuine quantum-information processing, since they all use 2 or 3 qubits. Yet the 
fact that all states — initial, intermediate, and final — that occur in a given NMR experiment 
are equivalent to ensembles of classical tops does not mean, by itself, that there is a classical 
model for the entire experiment. To see this, consider the following naive attempt to describe 
a unitary transformation U by classical transition probabilities: 

^p^\fi) =tT{UpU^Qi{n)) = j dn^wP{m)tT{UPi{Th)U^Qi{n)) . (3.9) 

= t[/(n|m) 

The transition function tu{h\m) is not a transition probability because it assumes negative 
values; for example, for the trivial case of the identity transformation, U = I, we have 

1 ^ 

tj{n\m) = tr(Pi(m)Qi(ri)) = 11(1 + 3m, ■ n,) . (3.10) 

The first criterion is the right one, however, for judging claims of having produced en- 
tangled states using bulk-ensemble NMR. For example, as pointed out in the claim of 
Laflamme et al. to have created a 3-qubit Greenberger-Horne-Zeilinger (GHZ) entangled 
state is incorrect. Despite the authors' assertion, "We describe the creation of a Greenberger- 
Horne-Zeilinger state of the form (|000) + |lll))/v^ (three maximally entangled quantum 
bits) using nuclear magnetic resonance. . . . We have thus extended the space of entangled 
quantum states explored systematically to three quantum bits ... ," no entanglement was 
created in that experiment; the statistics of any measurement performed on the purported 
GHZ state could have been explained in terms of classical correlations contained in w'^{n). 
Similar conclusions apply to the experiments in which Chuang et al. [|ll[] and Cory et al. |§] 
claim to have created 2-qubit Einstein-Podolsky-Rosen (EPR) entangled states. The states 
created in these experiment were unentangled, though Chuang et al. write, "As an applica- 
tion of the controlled-NOT gate, we used it in a simple quantum circuit to create entangled 
states from the thermal mixture. . . . We have experimentally confirmed this (nonclassi- 
cal) behavior, and the signature of the entanglement — a purely non-classical effect — is the 
strong reverse diagonal measured in the density matrix," and Cory et al. refer to applying 
the "spin-coherence XOR gate to a one-spin superposition to create an entangled state." 

Indeed, the first criterion can be applied to any NMR experiment that claims to have ma- 
nipulated quantum states in a particular way. For example, in describing the NMR version of 
quantum teleportation 0], a 3-qubit experiment, the authors assert, "Quantum-mechanical 
systems have information processing capabilities that are not possible with classical devices. 
One example is quantum teleportation, in which the quantum state of a system is trans- 
ported from one location to another without moving through the intervening space. . . . Here 
we report an experimental implementation of full quantum teleportation over inter-atomic 
distances using liquid-state nuclear magnetic resonance." This claim cannot be supported, 
because the quantum state at all stages of the experiment could be interpreted in terms of 
classical correlations among spin directions. What the experiment achieved was a reshuffling 
of these classical correlations. Nonetheless, as noted above, this conclusion does not imply 
that there is a classical model for the entire experiment. 
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IV. TRANSITION PROBABILITIES FOR QUANTUM GATES 



The essence of a classical model lies in the construction of nonnegative transition proba- 
bilities to describe unitary transformations. We first consider product unitaries, i.e., unitary 
transformations of the form V = Vi®- ■ - ^Vn, and then tackle the tougher task of entangling 
unitaries, i.e., unitary transformations that are not product unitaries. 

The expansion coefficients for the output state are given by 

w^P^\n) = tT(VpV^Qi{n)) = tT[pV^Qi{n)v) . (4.1) 

For a product unitary V, the unitary transform of Qi{n) factors into a product of transforma- 
tions for each qubit. If we introduce the 3-dimensional rotation operator Rj. (an orthogonal 
transformation) corresponding to the unitary operator \4, i.e., V^aVk = Rk(^, and use the 
fact that fik ■ Rk(T = Rk^fik ■ a, then the unitary transform of Qi{n) assumes the form 

V^Qiih)V = Qi{R^'nu • • • , Rn'^n) = QiiR-'n) • (4.2) 

Here R stands for the rotations on all the qubits. Now the output expansion coefficients 
assume the form 

yjVpV\n) = wP{R-^n) = w^iR^^ni, . . . ,i?^^n^) . (4.3) 

Not surprisingly, the action of the product unitary V is equivalent to rotating each of 
the classical tops. The corresponding nonnegative transition probabilities are w{n\m) = 
5{R~^n — m). This result shows that, for product unitaries, our transition probabilities are 
equivalent to the simple classical diagrams often used in the NMR literature . 

Before proceeding to entangling unitaries, we describe how decoherence is incorporated 
into our classical models. In some experiments decoherence is simulated as a mixture of 
product unitary transformations; an example is the use of gradient pulses in Naturally 
occurring decoherence processes act independently on the various qubits; on time scales 
short compared to the thermal relaxation time, they preserve the maximally mixed density 
operator. Any decoherence process that satisfies these two properties can be described as a 
mixture of product unitaries Thus decoherence can be handled in our classical models 
as a mixture of rotations of the qubits. 

To describe general, entangling unitaries, we need to introduce some additional notation. 
For an arbitrary density operator p and for < < 1, we define 

pe = {l- 9)M + dp . (4.4) 

Furthermore, we define the states 

Pe{n) = {l-e)M + ePi{n) (4.5) 

and associated operators 

i-e- ' 



Qe{n)^^—^I + e-'Q,{fi). (4.6) 
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For a general unitary operator f/, we define the transition probabilities 

wu{n\m) = ti {UP^{m)U^Qi{n)) = + rituin\m) > , (4.7) 

where tu{n\m) is defined in Eq. ( p.9|) . The nonnegativity of these transition probabilities 
follows from the argument leading to Eq. (|3.7|) , since U Pr,{fh)lJ'^ is a state of the form ( p.l|) 
with e = Tj. 

It is straightforward to write down the transition probabilities for elementary entangling 
gates. As an illustration, we compute them for the controUed-phase gate, 

U = C.,j = |0)(0|®1+|1)(1|®(T3 

= i(l + (73)^1 + ^(1-0-3)® (73, (4.8) 

acting on qubits i and j. We obtain 

wc,^{h\m) = ti (CijPr,{Th)CijQi{n) 

X ([1 + 3(mj)3(?2i)3][l + 3{mj)3{nj)3] + 9(mi x ni)^{mj x 71^)3 

+ 3[(mi)3 + 3(ni)3](mj_L ■ Uj^) + 3[(mj)3 + 3(71^)3] (mi_L ■ n^x) 

(4.9) 

where m_L is the projection of m into the 1-2 plane, obtained by setting to zero. As 
in this example, the qubits not involved in an entangling gate do appear in the transition 
probabilities, but in a simple, universal way. 

Applying the transition probabilities (|4.7| ) to an input ensemble = wP{n), we obtain 
a classical output ensemble 

Wontin) = I dVLrhWu{n\m)wiJffi) 



tr U 



n] 



Pv 



| + 77u;^^^'(n) . (4.10) 



(4vr) 

Unlike the quantum output ensemble the classical output ensemble suffers an 

increase in the fraction of molecules that are maximally mixed. This means that the mag- 
netization signal produced by the classical output ensemble is a factor of rj smaller than in 
the quantum description: 



J dQn'Wont{n){ni)f3^ ■ ■ ■ (niv)/3^ = tr(f/p.^f/V;3^ ® ■ ■ ■ ® ct/?^ 

= r]ti{UpU^ap^®---^af3^) . (4.11) 
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If we use the transition probabilities ( |4.7] ) to describe transitions of the spin directions in 
a gate-by-gate classical model of an NMR quantum computation, the magnetization signal 
of the model loses a factor of t] at each entangling gate. Though such a model gives a 
satisfactory account of an NMR experiment in which one ignores the strength of the signal, 
we can formulate a better model that avoids the decay for a few entangling gates. 



V. GATE-BY-GATE CLASSICAL MODEL 

The key idea in constructing an improved model is to introduce auxiliary, "hidden" spins 
ai, . . . , Oat = a, one for each qubit. We also need a "counter index" k, which increments by 
1 at each entangling gate. In the improved model, we can avoid the decay of the signal for 
K entangling gates, where K is the largest integer such that 

e < 77^+^ . (5.1) 

We assume that e < r^^ to ensure that K > 1. It is useful to introduce the function 

''' = [1, k>K. ('-2) 

Between the A;th and {k + l)th entangling gates (A; > 0), we represent a density operator 
p by the expansion coefficients 

= tr(pg,,^(a)) . (5.3) 



It is easy to show that these expansion coefficients are related to the original coefficients ( p.4|) 

by 

= + Vk'wTa) (5.4) 

and that p can be expanded as 

p = J dVL~aw1{a)P^^{h) . (5.5) 

For density operators of the form ( p.l|) , Eqs. (|5.4|) and ( p.7|) imply that 

1 _ e/„^'+i 

> -J^^ ^ (5.6) 

for all /c > 0. Throughout the experiment, the expansion coefficients w'^{a) can be interpreted 
as probability distributions. 

In our improved classical model, the gate operations manipulate the hidden spins; the 
effect of a gate shows up in the observable spins n through the correlations of n with d. 
Between the kth and {k + l)th entangling gates, this correlation is described by the joint 
distribution 

w^(n, a) = qkin\d)w^{d) , (5.7) 
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where 



qk{n\a) = 1^^^^ + r]kS{n - a) (5.8) 



is a conditional probability distribution. The marginal distribution for the observable spin 
variables, 

J dn-awl^h, ~a) = j^-p + Vkw'^in) = wf{n) , (5.9) 

is just the right mixture of w'^{n) with the uniform distribution to produce the distribution 
w'^ifi) that gives the measured expectation values. 

We now define transition probabilities for quantum gates, distinguishing as before be- 
tween product unitaries and entangling unitaries. For a product unitary V , we find, by an 
argument analogous to the previous one, that 

w]f^^\n, a) = w^(n, /2~^a) . (5.10) 

The effect of V is to rotate the hidden spins. Decoherence is handled, as previously, by 
mixtures of product unitaries. 

To deal with an entangling unitary transformation f/, we first note that the transition 
probabilities wu{a\b) defined in Eq. ( ^^ can be expressed as 

wu(a\h) = ti(UPr,e(h)U^Qe{a)) > . (5.11) 

Using this result, together with Eq. (|5l5|), we find 



p , 0<k<K 
Pn , k>K 

w^lf{h), 0<k<K, 
Wk+i (a) , k> K. 

To describe the effect of the {k + l)th entangling gate, we use transition probabilities 

Wk+i{n, a|m, h) = qk+i{n\d)wu{d\h) , (5.13) 
which take the input ensemble tUin = w1,{n, d) to an output ensemble 

.w(ft.a) = /<in,rfn,<«(n,a|A,i)u,„(A,6) 

= qk+i{n\d) / dVLiWu{d\h)wl(h) 



wS^r(n,a), Q<k<K, ^^^^^ 



w'llfinra). k>K 
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We have now constructed a gate-by-gate classical model in which there is no loss of 
magnetization signal for the first K entangling gates. How this works can be summarized 
as follows. The initial distribution Wo(o) of the hidden spin variables is chosen to have 
as little contribution from the uniform distribution as is consistent with nonnegativity. At 
each entangling gate, the ensemble of hidden spins suffers an increase in the proportion of 
molecules that are maximally mixed, but the observable spins retain the statistics of the 
quantum description by becoming more tightly correlated with the hidden spins. After 
K entangling gates, the observable and hidden spins become 5-correlated; thus from the 
[K -\- l)th entangling gate on, the signal decreases by a factor of rj at each entangling gate. 



VI. DISCUSSION 

Consider a bulk-ensemble NMR experiment whose parameter e satisfies the condition 
e < yf that underlies the argument in Sec. |V[ If we assume a = 2 x 10~^ in Eq. ( |2.2| ), this 
condition is fulfilled for < 6 qubits. The model constructed in Sec. provides an overall 
classical model in the sense of our second criterion: the unitary transformation U that maps 
the initial state to the final state is described in the model, with no loss of magnetization 
signal, by the transition probabilities (n, a|m, 6); mixtures of such transition probabilities 
are used to incorporate decoherence. The predictions of the model for the signal derived 
from the output state are identical to the quantum predictions. Since all NMR experiments 
to date involve 2 or 3 qubits, this overall classical model applies to all such experiments. 

We turn now to the imphcations of the gate-by-gate classical model constructed in Sec. 0. 
If we assume again that a = 2 x 10~^, in a 2-qubit experiment the model proceeds through 
K = b entangling gates with no signal loss, and in a 3-qubit experiment, through K = ?> 
entangling gates. We illustrate the implications by considering a particular 3-qubit exper- 
iment, the NMR version of quantum teleportation [Q. After preparation of a pseudopure 
state by the gradient-pulse technique pl, the teleportation experiment consisted of four 
operations: (i) an entangling 2-qubit gate, (ii) a 2-qubit mapping of the Bell basis to the 
computational basis, (iii) decoherence of two qubits in the computational basis, and (iv) a 
conditional 3-qubit unitary. Our model incorporates the decoherence step as a mixture of 
classical rotations of the hidden spins. Since the model can account for K = ?> entangling 
gates without loss of signal, it provides a gate-by-gate classical model whose predictions 
are identical to the quantum predictions. This experiment therefore does not demonstrate 
genuine quantum-information processing in the sense of our third criterion. For a larger 
value of a, when K = 2, our model predicts a drop in signal by a factor of 77 = 1/33 at the 
third entangling step. Nevertheless, even when K = 2, the classical model can still account 
for the teleportation experiment if the first two 2-qubit gates are "compiled" into a single 
3-qubit gate. 

Similar conclusions apply to the other 3-qubit experiments performed to date and 
to most of the 2-qubit experiments |jTl]-|l4|,||,|l5|,|l6| . The 2-qubit experiment described in 
ran through 15 entangling gates with an approximately exponential signal-to-noise decrease 
that is, however, much slower than the factor of r/ = 1/9 per entangling gate predicted by 
our model after the first K = 5 steps. The 2-qubit experiment reported in implemented 



up to 14 entangling gates, but the absence of signal-to-noise data makes comparison with 
our model difficult. 



12 



This paper, together with p3| , begins the task of estabhshing standards for assessing the 



quantumness of bulk-ensemble NMR quantum computation fS^JSJ] . Our gate- by-gate clas- 
sical model erects a hurdle in the way of NMR quantum- information processing. "Testing" 
our classical model is not the point, for no one would contend that it describes the physics 
that underlies an NMR experiment. The point is that experiments that fail to clear the 
hurdle can be explained in classical language and thus do no genuine quantum-information 
processing according to our third criterion. 

It will be easy for NMR experimenters to clear our hurdle, by measuring the signal-to- 
noise ratio in experiments involving many entangling gates, like the experiment reported in 
0. Yet jumping over our hurdle is probably not sufficient to provide a warranty of genuine 
quantum-information processing, for that would require showing that a given experiment 
is inconsistent with all gate-by-gate classical models, not just the model formulated in 
this paper. There are reasons for believing that our classical model is not optimal, the 
most cogent of which is that the model's transition probabilities for an entangling gate 
transform qubits that are unaffected by the gate. Avenues for improving the model include 
the following: (i) investigating more efficient representations both for quantum states near 
the maximally mixed density operator and for the transition probabilities between states, 
(ii) seeking gate representations that are more efficient when fewer qubits are involved in 
the gate, and (iii) addressing the extent to which one is allowed to "compile" successive 
entangling unitaries into a single operation (compiling is routine in NMR experiments as a 
method for reducing the length of a computation). 

We conjecture, however, that no matter how efficient the gate representations are made, 
an ultimate signal decrease is an unavoidable consequence of any attempt to describe en- 
tangling unitaries classically, even when the unitaries act only on separable states. Indeed, 
more interesting than our results would be a demonstration of this conjecture. Should the 
conjecture prove to be correct, one could conclude that the quantumness of NMR quantum 
computation, for small numbers of qubits, lies in the ability to avoid any signal decrease. 
More broadly, we speculate that the power of quantum-information processing comes not 
from entanglement itself, but rather from the information-processing capabilities of entan- 
gling unitaries p6|. 
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